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Abstract 

A field theory is developed for a tliermodynamical description of an array of par- 
allel non-singular screw dislocations in an elastic cylinder. The partition function 
of the system is considered in the functional integral form. Self-energy of the dis- 
location cores is chosen in the form suggested by the gauge-translational model 
of non-singular screw dislocation. It is shown that the system of the dislocations is 
equivalent to the two-dimensional Coulomb gas. The coupling potential is prevented 
from a short-distance divergency since the core energies are taken into account. Two- 
point correlation functions of the stress components are obtained. Renormalization 
of the shear modulus caused by the presence of the dislocations is studied in the 
approximation of non-interacting dislocation dipoles. It is demonstrated that the 
finite size of the dislocation cores results in a modification of the renormalization 
law. 
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1 Introduction 



Non-trivial topologically configurations in ordered states (e.g. vortices, dislocations, and 
other defects) attract appreciable attention in modern condensed matter physics. For 
instance, dislocations as imperfections of the crystalline ordering are of importance for 
structural, transport, and electronic properties of real solids. Two-dimensional ordered 
states are of special interest due to the significance of the topological defects for the 
corresponding phase transitions (THH]- The ideas of jlJ46j have been further elaborated 
for the description of the dislocation- mediated crystal melting in two dimensions [714TT]. 
The textbooks [T2|[T3] summarize an original approach to the ordered states and phase 
transitions dominated by the line-like disturbances. The peculiarity of [T2|IT3] is due to a 
systematical usage of (singular) gauge fields. Certain aspects of the statistical physics of 
dislocations are considered in [14r:18]. 

Dislocations have recently attracted considerable attention as far as the physics of the 
nanotubes is concerned |19fl22] . The point is that multilayer nanotubes can contain within 
their walls screw dislocations lying along the axis of the tube [21j. The external diameter 
of the nanotubes ranges from nanometers to tens of nanometers. The wall thickness 
varies from one to several tens of atomic layers. Therefore, effects influenced by the 
dislocation cores look attractive for study in the context of the nanotubes since the core 
sizes are comparable with those of atomic layers. Note that the electronic and mechanical 
properties of the graphene sheets in presence of dislocations are also of interest p3ti25] . 

According to the elasticity theory, the stress tensor components of a single dislocation 
are singular on the defect line. In reality, the stress components are smoothed out because 
of the formation of the core regions. Since the first attempts [2B|l?f]. various approaches to 
dislocations with non-trivial core are known. For instance, the quasi-continuum approach 
[28|l29]. the gradient elasticity [30,31j, and the Lagrangian translational gauging [331437] 
enable one to obtain appropriate solutions for non-singular screw dislocations. 

The elastic stresses of the screw dislocation are obtained in |33ti37j in the form of the 
superposition of the conventional far-reaching ("background") contribution and a short- 
ranged ("gauge") correction. The latter modifies the background fields within a compact 
core. Therefore, smoothing of the singular behaviour of the conventional screw dislocation 
occurs thus leading to the so-called modified screw dislocation. Note that the second- 
order elasticity is also of importance for the defects in crystals j2B]- The approach of [35j 
is suitable for study of the second order stresses in the case of a non-singular screw 
dislocation |39] . 

This paper is concerned with a thermodynamical description of a large number of the 
modified screw dislocations j35] lying within a long enough cylinder of a circular cross- 
section. The functional integration jlOl - HS] is used to represent the partition function and 
to study certain averages. The core energies are accounted for in accordance with [35j . 
It is shown that the collection of parallel non-singular dislocations is equivalent to the 
two-dimensional Coulomb-like system of charges interacting via the potential which is 
logarithmic at large distances but vanishes locally (smoothed out coupling). The two- 
dimensional Coulomb gas [iG] belongs to the class of systems covered by the theory [l -6j, 
and it demonstrates a relationship with the two-dimensional spin models [47ri49j . 

Dislocations available in a solid result in the renormalization of the corresponding 
elastic moduli [9l[lS|. Two-point correlation functions of the stress components are cal- 
culated in this paper for investigating of the renormalization of the shear modulus. It 



2 



is demonstrated that non-triviality of the cores results in a correction to the law of the 
renormalization of the shear modulus what could be appreciable for the nanotubes. The 
present field-theoretical approach is influenced technically by that of [SDIEI], where certain 
correlation functions have been calculated for the string models possessing the world-sheet 
vortices. 

The paper is organized as follows, section [1] is the introductory section. The partition 
function of the elastic cylinder containing the array of the modified screw dislocations is 
considered in section [2J Transformation of collection of the dislocations to a dual system 
of two-dimensional charges is presented in section |5J In section [31 we calculate certain 
thermodynamical averages, e.g., the average square of the dislocation dipole momentum 
and the stress-stress correlation function in the approximation of non-interacting disloca- 
tion dipoles. The renormalization of the shear modulus in the presence of non-singular 
dislocations is studied section |H Discussion in section [3] closes the paper. 



2 The partition function 

Consider a circular cylinder containing straight non-singular screw dislocations lying along 
the axis of the cylinder. The cylinder's material is approximated by elastically-isotropic 
continuum described by linear elasticity. Recall that the screw dislocation is characterized 
by the parallelism between its Burgers vector and tangent to the dislocation line |52[l53]. 
The modified screw dislocation j35j is a point of departure of the present investigation 
that treats the distribution of the dislocations as a thermodynamical ensemble at non- 
zero temperature. The functional integration approach is used below to investigate the 
partition function and certain correlation functions. In what follows, the Cartesian axis 
is along the cylinder's axis. 
Let us begin with the functional integral representation of the thermodynamical par- 
tition function Z of the elastic cylinder containing the non-singular dislocations: 

e-^'^V{a%,a\^,u,,e,,), (1) 
W = E- zEe^t , E = E,i + E,o,e , (2) 

where P is inverse of the absolute temperature T (the Boltzmann constant is unity). The 
functional W ([2]) consists of the following contributions (indices repeated imply summa- 
tion): 

Eel = + 

Ecore = j {leij{mce)ij - cr^^) d^x , (3) 

^ext = \j (y%{diUj + djUi-2Vij)d?x, 

where (inc e)ij denotes the double-curl operator: —eikiejmndk dmein i^iki is totally antisym- 
metric tensor). The functional E^i ([3]) is the elastic energy of superposition of two stresses, 
cr|^- and a^j, provided ^ and u are respectively identified as the shear modulus and the 
Poisson ratio. The notation a]', is reserved for the long-ranged contribution, while a^, de- 
scribes a non-conventional stress which modifies the background one, cXj , within the core 
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of the modified dislocation [35j 0. The dislocation core energy is given by -Ecore dS]) , where 
ejj(ince)jj originates (by means of linearizations) from the Hilbert-Einstein Lagrangian 
proposed in |35] in the framework of the translational gauge approach to dislocations H. 
Here Cij is the total strain tensor, and £ is a scale of the dislocation core energy. The term 
Eexi © is linear with respect to (symmetrized) derivatives of the displacement vector 
Ui, as well as with respect to a "source" Vij which is related to the plastic strain efj as 
follows: Vij = efj + Cij. The specific configuration of singular dislocation lines is pre- 
scribed by appropriately chosen since the plastic strain is concentrated on cut surfaces 
bounded by the dislocation lines [5^I5U]. The present approach enables one to avoid the 
stress divergencies on the dislocation lines. In this respect, an auxiliary field Cij (which 
is not a variable of the functional integration) is postulated to ensure the background 
fixing at ^ 7^ in such way that the components of afj just coincide with the conventional 
dislocation stresses governed by efj [Cij vanishes at £ = 0). 

The notation V{afj, a^^, Ui, Cij) implies the integration measure in ([T]) equal to a prod- 
uct of particular measures for each functional variable marked inside the brackets. More 
details on the definition of the functional integration measure by means of the appropriate 
limiting procedure can be found in [10] - H5] . The normalization factor is to absorb 
physically irrelevant multiplicative infinities. The partition function given by ([T]), ([2]) 
naturally generalizes the partition function of the array of singular dislocations [12ul3j . 
We shall apply the steepest descent approximation [lOl - HS] to Z given by ([1]), ([2j). The 
corresponding saddle point solutions are determined from requirement that the variation 
of the functional W vanishes, 5W = 0, under independent variations of its functional 
arguments. 

Firstly, the following observation concerning 2 ([T]) is in order provided three-dimensi- 
onal Eqs. (|3]) are viewed at £ = 0. Shifting subsequently a^j — )■ a^j + aij and Cij — )■ Cij + iij 
(here, aij and iij are new adjustable fields), one can remove in E the terms linear in a^j 
and Cij. Then, the integrations over a'^j and Cij are decoupled and can be compensated 
by redefinition of Eventually, E is reduced to E^i = Egila^'A, and two remaining 

integrations, over u and afj, result in the partition function of the array of singular 
dislocations characterized by a specific choice of [12[[T3]. 

Our framework is that of the plane elasticity since it is assumed that the cylinder is 
long enough |52 | l53 ] . The independence on the third coordinate reduces our study to the 
two dimensional problem. Therefore, the functionals ([3]) take the form: 



2yU _ 

^core = 2 j{iei{mce)i - a,') d^x , (4) 
Eext = / a^idiU-2Vi)d^x, 



where the integrands are (xi, 0:2 )-depen dent, and summation goes over i = 1,2. Since the 
displacement vector of the straight screw dislocation is along OX3 [52j, we use u = u^. In 



^ Two stress fields crj^- and a^^ correspond to splitting of the total elastic stress into, so-called, "back- 
ground" and "gauge" parts [34H37] . 

^The gauge Lagrangians quadratic in the dislocation densities are used in |341l36j what is equivalent, 
in the case of the screw dislocation, to the approach of [3S]. The most general three-dimensional gauge 
Lagrangian has been discussed in [SI] . Certain developments in the gauge gravity P^[551155] involve the 
gauge-translational features of the theory of dislocations. 
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addition, we abbreviate: a* = a* is b or c), Cj = e^s, etc. For instance, the notation 
for the functional integration measure is of the form: T>{af, crf, u, Cj). 

Let us consider equations (jl]) at £ = 0. The source Vi is taken at £ = in the form: 
Vi = ef . A single positively oriented straight screw dislocation intersecting the plane 
Xi0x2 at X = y (we use x = (xi,X2)) is "produced" by the plastic strain with a single 
non-zero component ef = {h/2)6{—xi + yi)S{y2 — X2) [59l|60]. Here, 9{-) is the Heaviside 
function and b is absolute value of the Burgers vector b lying along Ox^. Furthermore, 
the integration over u is equivalent to the insertion of the delta-functional 6{dia]^). The 
equilibrium equation dia]^ = is fulfilled by the Kroner ansatz = usij djf^ (see |61]), 
where eij is totally antisymmetric symbol of second rank. Vanishing of the variation 6W 
under the shift f^^f^+ Sf^ allows one to determine the potential 

Consider a collection of JV screw dislocations with the coordinates x = y/ and the 
Burgers vectors b/, 1 < / < A^. A disc of radius i? is a cross-section of the cylinder 
containing the dislocations. Assume that the defect positions on the disc are well separated 
from each other as well as from the boundary (the distribution of the dislocations is not 
too dense). For our purposes it is appropriate to consider the solution corresponding to 
the infinite plane, i.e., the influence of the boundary is neglected at a large enough R. 
Taking into account arbitrariness of 6f^, we obtain the extremum condition: 

M 

A/^ = 2^i^,e^ - S^ef ) = -zp(x) , p(x) ^ ^ j(x - Yi) , (5) 

7=1 

(2) (2) 

where j(x — y/) = 6/ 5 (x — y/), S (x) is the delta-function on XiOx2, and p(x) is thus 
the dislocation density. We determine from ([5]) using the Prandtl stress potential [6T] 
for the collection of the screw dislocations: 

= 7^J2^i^og\^-yj\. (6) 

7=1 

Now we turn to £ ^ and require vanishing of 5W for independent variations of 
the functional arguments of W. Generally, provided that Eqs. ([3]) are considered, the 
variation of af, results in the 'strain-stress' constitutive law: 

Also, the variation of Cij gives the equation: 

{mce)ij = ^a^j, (8) 

which is a linearized version of the Einstein-type gauge equation investigated in the non- 
linear approach of [_39j. Since di{mce)ij = 0, we put (in the plane problem) = usij djf^ 
to fulfil dia^ = 0. Arbitrariness of the variations and 6f^ results in the extremum 
conditions: 

A(/" + r) = 2i{diV2 

A(/' + r) = 



d2Vi) , 

2 _ 

n = ^ 



(9) 
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where Vi = ef + Cj, and Eqs. ©-([H]) are accounted for. Self-consistency of Eqs. ([9]) 
requires coincidence of their right-hand sides and takes the form: 

diC^ - d^C, = 1(p-zkV'=), (10) 
where p = p(x) is defined in ([5]) We represent a general solution to (fTOj) as follows: 

= '-{d,<P + e,,d,^), (11) 

where is the arbitrary regular function, while ip is to be found after substitution of (fTT]) 
to (fTOj) . If (fTOj) is respected, either or can be fixed arbitrarily. Our strategy ai i ^ 
is to keep still respecting Eq. Then Eqs. lead to a single equation for /^: 

(A - ^')r = ^p(x) , (12) 

which is solved as follows: 

_. AT 

= i^Y.^iKM^-yi\), (13) 
1=1 

where Kq{-) is the modified Bessel function. The choice 4> = 0, ip = f'^ ensures coincidence 
of (do]) and 1^. 

Therefore, ^■(/'^ + /^) is the stress potential of the collection of non-singular screw 
dislocations. The fields a]' and are rather the saddle point solutions of the functional 
W ([2]), while the proper stress fields are given by ^a^ and la^. Equations and (fT3|) 
demonstrate that a modification of the Prandtl stress potential ^f^ occurs within the 
cores, i.e., within the tubular vicinities of transverse size ~ k \ When k|x — y/| ^ 1, VJ, 
the total solution is dominated by the conventional contribution ^f^. The sum \{f^ + f^) 
behaves smoothly when x is close to either of y/, and the known singularities in the total 
stress distribution j^df + (jf) do not appear. 

Thus, the stationary solution corresponds to the modified screw dislocation [35] which 
is, in turn, in agreement with the solution obtained by means of the gradient elasticity 
131 ( 132 ] . The gradient elasticity itself belongs to a class of the generalized continuum 
theories which effectively take into account interatomic forces for the explanation of the 
material behaviour on the nano-scales (and thus inside the defect cores) [62j. The solution 
obtained will be used to estimate the partition function given by ([T]), and using 
the steepest descent approximation. 

Let us proceed with the consideration of the partition function Z. Integration by parts 
in W transforms ei(ince)j into the quadratic expression, —^{diCj — djCi)'^, invariant under 
the shift Cj — )■ Cj + dig, where g is arbitrary function (the Abelian gauge transformation). 
To restrict the gauge arbitrariness related to the functional integration over Cj, let us 
impose, so-called, "Coulomb gauge" diCi = [40j. Finally, the contribution £ej(ince)j in 
W is replaced by iciAet, and the shift 

e, ^ e, + {2£)-'A-'a^ (14) 

cancels the term linear in Cj. The Green function of two-dimensional Laplacian A^^ acts 
as convolution operator. The resulting Gaussian integration over e^. 
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is absorbed into and the partition function Z takes the form: 



Z = \- I e-^^V{alalu), (15) 



W = E - ^Eext , ^ = ^el + ^core , (16) 



where E^xt is given by (jlj), and -Ecore is expressed as follows: 

^core = ^/ ^-^-'^-d'x. (17) 

The functional W includes the elastic energy of non-singular dislocations Ed, while the 
core energies are given by (fT7|) . The representations expressed either by ([1]), ([2]), (jlj), or 
by (fT5|) . (fT6|) . (fT7|) are equivalent, i.e., the same Eqs. (j5]) and (Q arise for /"^ and f^, 
respectively. 

We are mainly interested in the exponential that contributes to the approximate ex- 
pression for Z (US]): Z ~ (DetA) e'^^^, where A is a kernel of the integral operator 
related to the second variation of W |13]. To derive W, it is necessary to substitute the 
solutions into W (fT6|) and apply the Green theorem. We neglect the terms expo- 

nentially small at kR ^ 1, as well as those contributions which are irrelevant provided the 
differences |y/j| = |y/ — yj| are moderate in comparison with R. The "electro-neutrality" 
condition bi = prevents divergency caused by a large log R. Finally, the effective 

energy W takes the form in the case of Af dislocations: 

W = 2fiJ 9ie^(x) ((A - K^y' - A"') diesis) d\ d^s , (18) 

where the plastic sources are expressed like in ([5]), i.e., (9ief = — p/2. Transforming 
W (ITS]) further, we represent it as the potential dependent on the dislocation positions 
{y/} = {y/}i</<A^: 

>V = >V({y/}) = 

-^/p(x)[log|x-s|+iro(.|x-s|)]p(s)d^.d^. 



b,h [log \J, - Jj\ + Mn\j, - jA)] + log(l-) Y. t?) 



where W(-) is defined as 



U{s) = log(|s) + Ko{s) , (20) 



and log(^) is the core energy of the dislocation in the Coulomb gas representation. The 
distances in (fT9|) are referred, for simplicity, to unit lattice spacing of a cubic crystal. 

_ dW 



We turn to calculation of the force F = — , . Its expression. 



-f = f{i^\yij\). fis) = ^ = ^ ^i(^)' (21) 



demonstrates that W (lT9|) at large separation |y7j| corresponds to the energy of the 
Coulomb attraction between two charges of unlike signs [bjbj < 0). The maximal attrac- 
tion of two opposite dislocations occurs at |y/ — yj| ?a 

Consider a thermodynamical ensemble of positive and negative modified screw dislo- 
cations located, respectively, at {y/}i</<Ar and {yj}i<i<j\f and possessing unit Burgers 
vectors. By analogy with the case of the two-dimensional electro-neutral plasma of pos- 
itive and negative charges [46j, the corresponding grand-canonical partition function Ze- 
is written, with regard at (I19p . as follows: 



+ (E ^(^\yi -yjD + H ^('^lyJ - yJl) - 2 w(«:|y/ - yji 



lytj i+j li^j 

(22) 

Here, 2M is the number of dislocations, and A is the chemical potential per dislocation. 
The representation (122|) implies that the summation over possible positions of the disloca- 
tions is replaced by the integration. The potential W(s) f l20|) is present in fl22|) instead of 
the regularised logarithm used in [50 » ,51 j for the case of the world-sheet vortices. It is ap- 
propriate to recall that the Coulomb gas of point charges is related to the two-dimensional 
sine-Gordon field theory 



3 The correlation functions 



3.1 Field-theoretical derivation of the stress-stress correlation 
function 

Gaining the experience of the derivation of the partition function (122|) . we pass on to the 
investigation of the relevant correlation functions. Let us take the field theory presented 
in section 2 as a starting point. Define the two-point stress- stress correlation functions 
{af{'Xi) (T*(x2))p (# is either b or c) by means of the following functional averages [40n45j : 

(af(xi)af(x,)), ^ i-y'af(x0af(x2)e-^^P(a^a^«,e,), (23) 

where the functional W is expressed by ([2]) and and the index P points out that the 
correlators (123|) are defined with respect to a prescribed dislocation distribution specified 
in terms of the plastic strain ef {i = 1,2). The latter is accounted for through the 
dependence of the functional W = W{Vi) on the source Vi. In turn, the normalization 
factor Zo is given by ([T]) though with W taken at = and i = 0. In other words, Z^, 
is the partiton function of a defectless state. 

Approach of the generating functional is a natural way to evaluate fl23|) in the frame- 
work of the functional integration. We introduce the generating functional ^[J^, J'^ | T*] 
as the functional integral parametrized by two auxiliary sources J*' and J'^, as well as by 
the source 
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Here the normalization factor is included into the integration measure, and each bold- 
faced notation, J* (7^ is b or c) or T', implies two components, say, jf and jf. The 
source is standing separately in left-hand side of f|24p since its meaning is different. 
The correlators we are interested in, (cr*(xi) cr*(x2))p, arise as follows: 

(af(xi)af(x2)), = Jim fc-^[J^J-|Pf] j^^f ^'^' g[J^J-|pP^]) , (25) 

jb.J'^^O V 5Jf (Xi)(5jj(x2) ^ 

where T'P^ is substituted for (= e-'^ + C) in order to express that Q is taken for a 
definite parametrization of the background configuration in terms of the field C (the 
appendix provides the calculation of ^[J*^, J"^ | T'^'^]). Moreover, T*"^ = 'P^^\^p_q-, i.e., the 



normalization is taken with respect of a defectless state, as in 

We shall obtain the correlation function (cr*°*(xi) cr*°*(x2))p of the physical stress field 
cr*°*(x) = 4(0"^ (x) + (Tj^(x)). First, we calculate the functional by shifts of the inte- 
gration variables which cancel the terms linear in J*^, J'^. It is appropriate to obtain the 
answer in terms of a single source Ji(x) = Jf{^) = </f (x). The result is of the form: 

g[3,3\ PP^] = ^[0, I ] , (26) 

where 



(27) 



j7i — t/j 26- 



(the Green functions A \ (A — k^) ^ imply the convolution operators). Then the physical 
correlation functions arise as follows: 

(.."(xO^r (xd), = I- (e- [J. J I Tf] Jj£^^6lJ, J I P*]) . (28) 

We calculate (EH]) using ([26]) and (127]), and obtain: 

(^^(xOaf (X2))^ = ^e-'^^ lim[(5,^(.,)g)(5,^.(.,)g)--^5,^(.,)5,^,(.,)Q] , (29) 

where 5ji(x) = S/6Ji{'x), and W, Q are given by (fT9]) and (127|) . Equation ( 129]) is re- 
expressed after the variational differentiations: 

(ar*(xi)af (X2)), = e-/^^ k°*(xOaf (x^) - /-a(.,).a(.,)^.W(/c|xi - x^l)] , (30) 



27r/3 



where 



^Cifc 5(x)fe y W(k|x - s|) 9ief (s)ds (31) 

is the total elastic stress of the dislocational configuration. The exponential factor in 
( 150]) is the Boltzmann weight containing in its exponent the energy of the array of the 
dislocations. 

Consider, for a comparison, our theory without the influence of the core energy, i.e., 
at £ = 0. Derive, firstly, the stress-stress correlator {af{'Xi)a^{-X2)) in the absence of 



the defects. The source ef being treated as unphysical one (that becomes zero after the 
variation) enables us to obtain: 

(a,^(x0a^^(x2)) = ^ 9(.,).^(x.). log |xi - x^l . (32) 

Let the plastic source be the physical one. Then we obtain: 



(a,^(xi)or)'(x2))p = e ^(9(xi),5(x2)^. log|xi -X2I 

+ ^ XI blbj{^iA^ih log 1^1 - y^l) {^jl9{^2)i log |X2 - yj|) 
I,J 



(33) 



where W is given by (|T9|) . although with U replaced appropriately by the logarithmic 
potential. 

The representation (150]) corresponds to a specific spatial distribution of collection of 
the dislocations. We shall average (cr*°*(xi)cr*°*(x2))p over positions of the defects. 

3.2 Mean square of the dipole momentum 

Let us investigate the grand-canonical ensemble of the dislocations using the so-called 
dipole phase approximation which corresponds to pairs of the dislocations with opposite 
signs (located in xf, 1 < I < Af) bound into "molecules". The corresponding partition 
function is specialized as follows |50] : 



00 ^ M 



M=Q 1=1 ■ 



( 



where w{r}j) is the energy of dipole centered in = {x.f + x7)/2 with the dipole 
momentum rjj = x.'j' — x.J , while wjj is the energy of interaction between J*^ and J^^ 
dipoles. 

Firstly, we calculate the mean square of the dipole momentum for a single molecule: 



exp{—l3w{r])) r] dr] 

{V') = ■ (35) 

exp{—f3w{ri)) drj 



The average f l35|) has been calculated in |3] for electrically-neutral gas of particles with 
charges ±g interacting through the potential: 

Ixj — Xjl 

2A - 2qj qj log , x/ - xj > a , 

a 

0, |x/ -xjI < a, 

where x/ is the charge position, and a is the cutoff (e.g., the particle diameter, or 
lattice spacing). Moreover, 2A is the energy necessary to create a pair of opposite charges 
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at the distance a. The dipole energy is given by Pw{r}) = K. \og{\ri\ / a) , where K, = 2/3g^ 
The average ( 15^ takes the form |3]: 



provided that /3g^ > 2. Therefore, (ry^) ^ c? in the hmit of zero temperature, while (77^) 
infinitely grows at the temperature T near its critical value given by fic'f' = 2. 

Two-particle potential studied in [50j is expressed as (3w{r]) = 2tt}CU{\t]\), where 
/C = 2/3g^ (here /3 is an effective parameter viewed as inverse temperature), and V({ri) is 
regularised at small t] = \r}\, Uirj) = | log " . We obtain from (135|) at /C > 2/7r: 

(37) 

The dipole momentum (1371) is strictly smaller than (77^) (15^ at /C > 4, i.e., in the second 
case the molecules are more compact at small enough temperatures. 

The present paper deals with the dipole energy U (120|) . and the average f l35|) is specified 
as follows: 

exp(-/C(logr7 + KQ{r]))) rf dr] 

^^\rn=-roo , (38) 

exp(-/C(logr7 + Ko{r]))) 7] dr] 



c 





where /C = ^}P'P/2tx. The integral in the nominator of (138|) diverges at /C < 4. Since 
dissociation is most probable for pairs of the dislocations with unit Burgers vectors, the 
dipolar phase does not exist at the temperature T > = 

It has been proposed in |31 ] IM | l36] to represent the radius of the dislocation core r, 
in the form — Tjb/K, where rjb = 4.0 [3T], rjb = 6.0 [36], or rjb = 10.0 j3l]. In turn, 
K has been selected in terms of interatomic spacing a as follows: k, = 0.25 a [31j, 
K = 0.399 a [SniES]- For instance, Tc ~ a according to |31j . while ~ 2.4a according 
to [SniES]- We represent the ratio ( 155]) at /C close enough to its "critical" value as follows: 

2r)^ 1 -1 P'' 

- -TT^^n-r: T—. ^ A" = / exp(-4(logr/ + iro(r/)))r/dr/, (39) 



ft:2(l + 2r/2A"') /C-4' 

where rjb and k are taken according to [30], [3T] or [36]. Moreover, double-sided estimates 
can be obtained for tl{ri) ( 120|) by adjusting appropriate trial functions. Such estimates 
enable one to demonstrate that dependence of {rf) (155]) on growing K, is characterized by 
the double inequality: 

h^"' (^') < ^ • (40) 

The potential Uijf) = | log(l -|- r]"^) can be used in (138]) instead of U ( 120]) for a comparison 
of extensions of the dipole momenta. According to the estimate (140]) . the dipoles are more 
compact for U ( 120]) at small enough temperatures. 

3.3 The correlator of two stresses 

Now, we turn to the averaging of the correlation functions ( 130]) over the dipole positions. 
We shall follow [50j where the asymptotic behaviour of the correlation functions of the 
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string models have been investigated in the presence of the vortex dipoles. We shall denote 
the new average as ((cr*°*(xi) 0"*°*(x2))), and its expression arises as follows: 

K°*(xi)af (X2))) = ^5(.o.^{x.).W(«:|xi-X2|) 

pw (41) 



numbers of dipoles, 
dipolo positions 



where Zdip is the partition function flM|) . The dipole positions are confined within a disc 
of radius R. 

According to f HO]) . the dipoles are very compact since the dipole momenta are not 
too large at small enough temperature: {rf) <^ /t^^. Therefore, summation over the 
dipole positions can be replaced by integration. We use the dipole's center of mass and 
momentum coordinates, respectively, = (y^ + y]^)/2 and 77^ = yj — y^ for each pair 
of opposite dislocations at y^ forming the V"^ dipole {1 < L < TV). We continue to use 
the approximation of non-interacting dipoles, and the sum in right-hand side of ( HTj) takes 
the form: 



E ar(xi)ar(x2)e-^^ = 

numbers , positions 

X f (W(/s:|xi - y+ 1) - W(«:|xi - y^|)) {U{k\x, - yt\) - U{k\^2 - I)) 

K,L=1 



(42) 



where M is the number of dipoles, and we use for cr*°*(x) the representation fl3Tl) . We 
estimate \rii\ <^ |x — (compactness of the dipoles) and adopt in leading approximation 



W(K|x-y+|) - W(K|x-y^|) ^ -(r7^,9jW(«:|x-|J), (43) 

where (-, ■) stands for the scalar product of 2-vectors. The relation ( 1431) allows us to 
re-express fl42p as follows: 



/3(2A + «'(r?^)) 



(44) 



M=l ' 1=1 

M 

X E (^L, <9xJ W(/«|Xi - ^^1) (77^, (9xJ W(k|x2 - ii}) . 
L=l 

To proceed with essential technical task is to calculate the integral: 

j I d2r7(r7,9.jW(«:|xi-||)(r7,9xjW(«:|x2-^|)e-'^(2^ + '"(^)). (45) 



First, we express the ry-integration by means of the relation 

exp(-2/3A - /C(log(|/«r7) + Ko{Kr])))r]ir]j d^ij 



7 

Hri^). (46) 



j exp(-2/3A - /C(log(|«:r/) + K,{kt^))) d\ ^ 
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where /C = 1] = \r]\, and (ry^) is given by (138|) . Then, after introducing the notation 
N for the dipole density jHED]: 

N = y"exp(-2/3A-/C(log(|«:r/) + iroM))d2r7, (47) 



the relation (146|) allows us to re-express the integral (145|) as follows: 



y W(k|xi - id , W(/«|X2 - ID) d^l , 



(48) 



The Green theorem enables one to carry out the ^-integration in ( HHj) . Eventually, we use 
the partition function Zjip and obtain: 

((ar(xO<Hx2))) = ^9(.,).a(.,)^W(«:|AxD 

- u^fo^ kIAxI 
- {ri')N ^ e,ke,i ^(^^.^{x^), (W(«:| AxD A^i(/t| AxD) , 

where (and henceforth) Ax = xi — X2. 

At large separation of the arguments, |Ax| 3> n , the asymptotic of fl49p is governed 
by the logarithmic contribution since the modified Bessel functions decay exponentially. 
We use the rule of differentiation of the logarithm, 

5(xih<9(x2hl0g|Xi -X2I = 

1 / , JXi-X2)fc(Xl-X2)z ^ (2) (50) 

-cifei + 2 . ^5 - ndki d (xi - X2) , 



|Xi - X2I" \ |Xi - X2 

as well as the relation ejfcCji = Sij6ki — SuSkj, and obtain: 



2(Ax), (Ax),- 




Ax 


2 




Ax 


4 ' 



((ar(xi)a-Hx2))) = [^ + ^(^^) 

Therefore, the stress-stress correlator (being considered with dislocations as well as with- 
out them) decreases as |Ax|~^ at growing separation |Ax|. 

The asymptotic law |Ax|^^ is due to the analytical structure of the logarithm dif- 
ferentiated. This law being extrapolated to small distances, |Ax| ^ 1, would imply a 
singular behaviour of the stress-stress correlator in question. In the present approach, 
the short-distance behaviour of the correlation function is less singular because of the 
core influence. Since the logarithmic contributions are canceled in the potential U, the 
short-distance behaviour is also less singular: 

((a- (xO o--(x2))) = ^ Y [- {I + N{ri') ^ - log(|/.| Ax| 

(52) 

+ ^^^|^ + 0(|Axriog|AxD]. 

We considered the approximation of non-interacting dipoles. Dipole-dipole interaction 
can also be taken into account in the present framework. Moreover, second order correc- 
tions to the stress fields inside the core could become important for sufficiently dense gas 
of the dipoles [39]. However, this should be a subject of a separate investigation. 
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4 The renormalization of the shear modulus 



The present section is to investigate the renormalization of the shear modulus caused 
by the modified screw dislocations. One should refer to the original paper [9j for the 
definition of the inverse of the tensor of renormalized elastic constants in the presence of 
dislocations. Further details on the dislocation contribution to the elastic constants can 
be found in [5] and [TB] (devoted, respectively, to two- and three-dimensional situations). 
Adopting the relevant definitions jSlIIS], we shall consider the following expression for the 
renormalized shear modulus fj,ren- 

^ = ^ E // ^"'(^2))) d^xid^x^ , (53) 

where the correlator is defined by (141 p . S is the area of the sample's cross- sect ion. We 
use (H^ and f l50p . and obtain: 

E ((^r(^i)^r(^2))) = ^^(^o(«:|Ax|) + 

fc=i.2 _ (54) 

+ ^k\Ax\Ki{k\Ax\)) , a= 



2 ' ' ' 'V 2 

"Non-diagonal" correlators (((T^°'^(xi) cr*°*(x2))), k ^ I, are negligible with respect to the 
two integrations in (153|) . The parameter a (15^ is proportional to the mean area covered 
by the dipoles, 2TT{rf)N . Therefore, one obtains from ( 153|) and ( 154|) the following answer: 

= -Ci{kR) + aC2{KR), (55) 

/^ren 

where the functions Ci{nR) and C2{nR) are given by the modified Bessel functions: 
Ci{kR) = I - 2Ki{kR)Ii{kR) , 

(56) 

C2{kR) = 2 - 2h{KR){Ki{KR) - kRK[{kR)) , 

and K[(z) = j^Ki{z). The renormalization rule fl55|) demonstrates the dependence of 
the shear modulus /iren on the dimensionless parameter kR. The coefficients Ci{kR) and 
C2{kR) both are positive and less than unity though tend to unity at nR — )■ 00. We obtain 
the estimates for Ci{kR), C2{k,R) at increasing kR: 

Ci{kR) ^ 1 - + . . . , C2{kR) ^ 1 - + . . . , (57) 

where the ellipsis imply the terms 0{{k,R)~'^). 

On the other hand, an analogue of (I54p valid for singular screw dislocations takes the 
form: 

E ((^'(^1) ^'(^2))) = ^(1 + «/i) ?(xi - X2) . (58) 

k=l,2 ^ 

The right-hand sides of (15^ and f l58|) are weakly coinciding at k ^ — (i.e., when 
the core's scale is shrunk). This can be demonstrated by means of integration with an 
appropriate trial function. Inserting fl58l) into fl53l) one obtains: 

= - + a. (59) 
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The renormalization rule (159|) is in agreement with that obtained in the original paper |1| 
with the help of the macroscopic stress function of collection of the dislocation dipoles. 
Equation (15^ agrees with the renormalization of the shear modulus found in [HlIIS] pro- 
vided the concentration of the defects is low. As is seen from (159|) . increasing of a results 
in decreasing of /iren- Clearly, Eq. (155|) is reduced to (l59l) when nR tends to infinity, and 
so the unit values of Ci and C2 correspond to the case of singular dislocations. Roughly 
speaking, shrinking up the core regions one goes back to singular dislocations. 
Equation (155]) can be re-expressed as follows: 



At finite nR (the scale 1/k is comparable with the sample's scale i?), the non-triviality of 
the dislocation cores is valuable for the renormalization of the shear modulus. According to 
(lFr|) and (I60p . the character of the decreasing of /Xren is changed in comparison with the case 
of the singular dislocations: it is slower since C2{k,R) / Ci{kR) < 1 at large kR. Recall that 
the critical exponents of appropriate correlation functions are analogously renormalized 
because of the presence of the vortex pairs, for instance, in the two-dimensional Bose 
gas [6j, in the classical planar Heisenberg model [47J, or on the string world-sheets |50 ^ i51 j . 

5 Discussion 

An array of parallel singular screw dislocations is equivalent (as a planar system with 
respect of the sample's cross-section) to the two-dimensional Coulomb gas of charged point 
particles. The latter is a subject covered by the theory [THT^ of the phase transitions 
in the low-dimensional systems of condensed matter physics. Dislocations influence re- 
normalization of the elastic constants of the corresponding material sample. This paper 
investigates the re-normalization of the shear modulus in the case of the screw dislocations 
possessing the core regions of finite-size. The transformation to the Coulomb-like system 
is used. 

The approach of |35 [ l39 ] to description of singularityless dislocations is elaborated in 
the given paper further for studying the collection of the modified screw dislocations as 
a thermodynamic ensemble. Specifically, a long enough cylinder pierced by non-singular 
screw dislocations is studied. A field-theoretical formalism is developed for investigating 
the corresponding partition function in the form of the functional integral. The stress 
potentials of the modified dislocations play the role of its saddle points because of the 
choice of the energy functional. The plastic external source is involved which governs 
the background stress distribution. This is distinct from |3lH37] where the usage of the 
so-called "nuU-Lagrangian" is responsible for arising of the background stress field. 

Calculation of the partition function relates the system of dislocation dipoles to equiva- 
lent description of the electrically neutral Coulomb gas of charges interacting via potential 
which is logarithmic at large separation but tends to zero for the charges sufficiently close 
to each other. The smoothing of the Coulomb potential at short mutual separations 
occurs since the self-energy of the cores is accounted for. The stress-stress correlation 
functions are obtained and used for studying the renormalization of the shear modulus in 
the approximation of the dilute gas of the dislocation dipoles. 

It is demonstrated that the renormalized shear modulus depends non-convent ionally 
on the ratio R/ k = kR of two lengths characterizing the sample's cross-section and the 




(60) 
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dislocation core sizes. Note that applicability of the effects of renormalization of the elastic 
constants to experimental observations is discussed in [TB] for the case of singular defects. 
Since the contributions caused by the core are sensible at moderate kR, it is hopeful that 
the formalism developed could be efficient for the nanotubes with comparable R and k . 
With regard to [32l[35l|37] , it is hopeful that the formalism presented can be extended to 
a hollow cylinder, as well as developed further to cover the modified edge dislocations. 
The latter could be interesting as far as the physics of multi-layer nanotubes and wrapped 
crystals is concerned 
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Appendix 



Consider the generating functional ^[J^jJ*^!?^] ([Mj) dependent on three 2-component 

b 



sources, J , J^, and P: 



The exponent in (Al) is specified, after fixing the "Coulomb gauge", as follows: 
-PW + 1(3 j {J^a]' + r^a]) d?x = 



2fi 



{a'^ + a^Yd^x - 2(3i I e.Ae^d 



X 



+ zf3j a]' {diU - 2Vi + J^) d^x + (3 J {2ei + iJ^) d^ 



X . 



It is appropriate to calculate (Al) by shifts of the functional integration variables 
As a first step, the strain field Cj should be integrated out by the shift 



Ci + 



21A 



After re-arrangements we obtain for ^[J*^, J*^ | T']: 



^[J^, J"|P] = const X / exp 



zl 

2/i 



al D-^al d^x + i(3 / al d'^x 



(a^a^ - z2^,a^{d,u - 2V. + J^' + -al)) d^ 



X 



(Al) 



{A3) 



2fi J \ " " • " X ' 

Decoupled integrations, being constant factors which are not of interest now, are system- 
atically included into the prefactor const. The kernel is defined in (A3) as follows: 

{AA) 
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where A is the Green function of two-dimensional Laplacian, and 5 is the delta-function. 
Next step is to shift subsequently the variables and a'[. It is assumed that j}" = 
= Ji- After the shifts 

(^"^ (^^^ + (diU - 2ef + 
— cr'^ — ifi D{diU — 2e. 

the generating functional takes the form 



) 5 



Q[3, J I P] = const X / exp 



-2{d,u + J,){DJ, + 2^)) d'x -tl3 j g\{2C, - D(9,n - 2ef )) d^ 



X 



|- / (a,^ D- V,^ + (T,Vf + 2(T,Vn d"x 



(A5) 



(A6) 



where the operator D is inverse to (A4), and Ji = J,- 2ef . 

The experience of section 2 tells us about a necessity of fixing of the background 
contribution by means of a self-consistent choice of the functions Cj. From a viewpoint 
of the functional integration, the corresponding substitute 



a = - D {d,u - 2el 



(which solves (fTOj) ) allows one to decouple the integration over u. The corresponding 
functional integral is just responsible for the dependence of the partition function on the 
plastic strain ef and thus on the defect distribution. Then the following representation 
arises: 

g\i, J I PP"^] = const X j exp ^ j {{OiU + Ji){5 - D){diU + Ji) + J^DJ^) d\ V{u) 

m 

where the notation 7^^^ implies that the arbitrariness due to Ci is removed. As a final 
step, the shift 



u 



1 

A 



u -diJi, 



allows one to get rid of the contribution of the first order in diU. Finally, the generating 
functional takes the form: 



g[3,3\V^''] = ^;[0,0|Pf]x exp 



-Hl3 



A-k2 



Ji (fx 



A A-/t2 



dkJ^k 



where J,i = J,i — 2ef. All the integrations decoupled are gathered within ^[0,0 17^^^] 
where V"^^ = 'P'^^\^p_q- Equation (A9) is just the answer expressed by (12^ and (1^ . 
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